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We examine how light neutrinos coming from distant active galactic nuclei (AGN) and similar 
high energy sources may be used as tools to probe non-standard physics. In particular we discuss 
how studying the energy spectra of each neutrino flavour coming from such distant sources and 
their distortion relative to each other may serve as pointers to exotic physics such as neutrino decay, 
Lorentz symmetry violation, pseudo-Dirac effects, CP and CPT violation and quantum decoherence. 
This allows us to probe hitherto unexplored ranges of parameters for the above cases, for example 
lifetimes in the range 10~ 3 -10 4 s/eVfor the case of neutrino decay. We show that standard neutrino 
oscillations ensure that the different flavours arrive at the earth with similar shapes even if their 
flavour spectra at source may differ strongly in both shape and magnitude. As a result, observed 
differences between the spectra of various flavours at the detector would be signatures of non- 
standard physics altering neutrino fluxes during propagation rather than those arising during their 
production at source. Since detection of ultra-high energy (UHE) neutrinos is perhaps imminent, it 
is possible that such differences in spectral shapes will be tested in neutrino detectors in the near 
future. To that end, using the IceCube detector as an example, we show how our results translate 
to observable shower and muon-track event rates. 

PACS numbers: 14.60.Pq,14.60.Lm,14.60.St,13.15.+g,11.30.Cp,98.54.Cm,95.85.Ry 



I. INTRODUCTION 



The extraordinary success of neutrino experiments in 
the last few decades has propelled neutrino physics to 
the centre-stage of particle physics. A series of seminal 
observations [TlH2| have provided us with a "new stan- 
dard model" (nuSM), in which the standard model of 
elementary particles is augmented by three massive neu- 
trinos which mix. Therefore, in addition to the standard 
model parameters, the nuSM also includes at least 2 mass 
squared differences 1 Am^ and Ato^, three mixing an- 
gles #i2, #23 and #13, one (so-called) Dirac CP phase, 
and two Majorana CP phases (if neutrinos are Majo- 
rana particles). With the existence of neutrino masses 
and mixing confirmed, focus has now shifted to the next 
level, viz., (i) making precise measurements of the known 
oscillation parameters, and (ii) determining the hitherto 
unknown properties of neutrinos. Among the neutrino 
parameters which belong to the standard nuSM picture 
and which are still unknown are the mixing angle #13, 
the sign of Amfi, and the CP phase(s). Next-generation 
neutrino experiments are expected to throw light on some 
or all of these standard neutrino oscillation parameters. 
Data from these experiments can also be used to probe 
physics beyond the standard paradigm. This includes 
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a variety of new physics scenarios such as non-standard 
interactions at the source and detector as well as dur- 
ing the propagation of neutrinos |13H15j . non-unitarity 
of the neutrino mixing matrix [16 18 , violation of the 
equivalence principle [HO [2U] , neutrino decay [2"TH2"5] , vi- 
olation of Lorentz invariance [2H US] , pseudo-Dirac neu- 
trinos [2"6Tt28| . and neutrino decoherence [25]. Among 
the forthcoming sets of experiments which hold promise 
for new physics searches are neutrino telescopes, which 
have been designed to observe ultra high energy neutri- 
nos coming from astrophysical sources. 

Very high energy cosmic rays with energies as high as 
10 10 GeV have been observed [301 EI]. There is also now 
a large body of evidence for high energy gamma rays 
coming from astrophysical sources. Understanding the 
origin and source of these high energy cosmic rays re- 
mains a challenge. A plethora of papers have appeared 
in the past trying to provide a viable model for these 
observations. Nearly all such models allow acceleration 
of protons to energies in the realm of 10 10 — 10 11 GeV. 
Such high energy protons will invariably lead to the pro- 
duction of highly accelerated pions (and kaons) through 
P7 and pp collisions. These pions would in turn pro- 
duce neutrinos carrying energy anywhere in the range of 
10 4 — 10 10 GeV depending on the type of source. Detec- 
tors such as AMANDA [32], IceCube [33J, BAIKAL [33], 
ANTARES [35], KM3NET [35], RICE [37] and ANITA 
[38) have been constructed (or are under construction) 
using techniques that would make it possible for them to 
observe these ultra high energy neutrinos. 

Neutrinos produced via decay of pions are expected to 
roughly carry the flavor ratio {v e : : v T =) 1 : 2 : 
at the source. Standard neutrino oscillations in vacuum 
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massage this ratio during propagation to 1 : 1 : 1 [SSI HO] 
at the detector, if we assume 8 13 = and 9 2 3 = 7r/4 
consistent with current data [4"IH4"3"] . It has also been re- 
cently stressed [H] that standard flavor oscillations over 
Mega-parsec distances make the neutrino spectra of every 
flavor nearly identical in shape. Therefore, if for any rea- 
son the astrophysics in the source leads to a ratio different 
from 1 : 2 : or spectral shapes for flavours which differ 
widely from each other, standard oscillations still mas- 
sage them into identical shapes and magnitudes which 
are within a factor of roughly 2 of each other by the time 
they reach the earth [33] . 

The potential of the neutrino telescopes to probe new 
physics using absolute flux ratios has been studied in 
[221 [45 - 53 . In an earlier paper [44] , we considered two 
specific new physics scenarios, viz., neutrino decay and 
Lorentz invariance violation, and showed how they affect 
the diffuse ultra high energy neutrino flux. We empha- 
sised how spectral information could be used to extract 
new physics from the ultra high energy neutrino data. 
In this paper we extend our earlier analysis to include 
more new physics cases. In addition to neutrino decay 
and Lorentz invariance violation we consider the effect of 
pseudo-Dirac neutrinos and neutrino decoherence during 
propagation. Further, we include three flavour effects, 
by allowing the mixing angle sin 2 (#13) to vary from to 
0.1 and the CP phase — 2n. Thus, we also take into 
account the uncertainties in the present values of these 
poorly known parameters. We will show that the uncer- 
tainties in these poorly known parameters cannot mask 
the effects due to neutrino decay or Lorentz-violation. 

We calculate the diffuse ultra high energy neutrino 
fluxes and for specificity focus on active galactic nuclei 
(AGN) as sources for this flux. We demonstrate how 
the diffuse flux spectra change as a result of new physics 
scenarios. Neutrino decay, depending of the choice of 
the neutrino lifetime, results in partial-to-complete disap- 
pearance of the heavier neutrino mass eigenstates, leav- 
ing mainly the lightest mass eigenstate to be recorded in 
the detector 2 . Therefore for complete decay, the flavor 
ratios at the detector are given by \U ei \ 2 : lU^] 2 : \U Ti \ 2 , 
(i = 1 or 3) which for tribimaximal mixing is 4 : 1 : 1 for 
the normal hierarchy (i — 1), and 0:1:1 for the inverted 
hierarchy (i = 3). However, it is possible that neutrino 
lifetimes are such that the decay is not complete and only 
occurs for the lower energy neutrinos. This would intro- 
duce an energy dependence which will change not just 
the ratios but also the spectral shapes. 

On the other hand, the effect of Lorentz invariance 
violation is more pronounced for higher energy neutri- 
nos. In particular, at higher energies a breakdown of 
Lorentz symmetry will lead to the breaking of the ex- 
act/approximate [i — t symmetry that exists for stan- 
dard neutrinos. In fact, we will show that for values of 
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the Lorentz invariance breaking parameter a > 10~ 26 
GeV, there are almost no r-neutrinos arriving at the de- 
tector above E > 10 5 GeV, whereas the !/„ flux is en- 
hanced compared to its expected values. Observation of 
this large breaking of the fJ> — t symmetry at higher en- 
ergies by ANITA or Auger would then be an indication 
of a possible breaking of Lorentz invariance. 

In this work we also estimate the number of muon 
track and shower events in IceCube to demonstrate how 
this method involving spectral distortions can actually 
be used in the terrestrial neutrino telescopes. We show 
the flavor ratios not just in terms of the diffuse fluxes but 
also in terms of the ratio of muon track to shower events. 
This event ratio is seen to have spectral distortions at 
the low energy end for the case of neutrino decay. 

The paper is organized as follows. In Section [TT] we 
briefly outline the procedure for calculating the diffuse 
UHE neutrino fluxes from AGN allowing standard oscil- 
lation among the flavours during the propagation of these 
neutrinos. In the next section, we demonstrate the effect 
of standard oscillations in even-ing out the shapes and 
magnitudes of the flavour fluxes from AGN. Section |IV| 
then shows the modification of these fluxes due to decay 
of the heavier neutrinos, and its effect on the number 
of detectable events at a large volume detector like the 
IceCube. We examine the effect of variation of #13 and 
the CP violating phase Sep in Section [V] We look at the 
effect of Lorentz-symmetry violation in Section |VI[ and 
finish with brief investigations of the effects of pseudo- 
Dirac neutrinos and decoherence in the last two sections. 



II. THE DIFFUSE NEUTRINO FLUX FROM 
ACTIVE GALACTIC NUCLEI 

Active galactic nuclei are extremely distant galactic 
cores having very high densities and temperatures. Due 
to the high temperatures and the presence of strong 
electromagnetic fields, AGN's act as accelerators of fun- 
damental particles, driving them to ultra-high energies 
(> 1000 GeV). The acceleration of electrons as well as 
protons (or ions) by strong magnetic fields in cosmic 
accelerators like AGN's leads to neutrino production. 
Specifically, accelerated electrons lose their energy via 
synchrotron radiation in the magnetic field leading to 
emission of photons that act as targets for the accel- 
erated protons to undergo photo-hadronic interactions. 
This leads to the production of mesons which are un- 
stable and decay. In the standard case the charged 
pions decay primarily contributing to neutrino produc- 
tion via tt^ 1 — > fi and subsequent muon decay via 
\ir — > e VpV,,. This leads to a flavour flux ratio of 
(y e : : v T =) 1 : 2 : in the standard case. The parti- 
cles finally produced as a result of this process are, thus, 
high energy neutrons, photons, electron pairs and neu- 
trinos. 

In this section we calculate the diffuse flux spectrum 
of neutrinos escaping from both optically thick AGN's, 



3 



which are so called because they are opaque to neutrons 
and trap them, and optically thin AGN's, which are 
neutron-transparent, and detected at distant detectors, 
for instance, at IceCube [32] . 

To calculate the flux from optically thick sources, we 
use the spectra of neutrinos produced in a standard AGN 
source, as discussed in detail in We then account 

for red-shifting in the energy dependence of the spectra 
appropriately. To obtain the upper bound for the dif- 
fuse AGN flux we vary the break energy Eb within the 
allowed range and maximally superpose all the result- 
ing spectra. To obtain the diffuse AGN flux spectrum 
at earth using a standard AGN distribution across the 
universe, we integrate the red-shifted spectra from the 
individual sources over the standard AGN distribution 
in the universe. The resulting diffuse bound and spec- 
trum are then normalised using the cosmic ray bounds 
also obtainable using a similar calculation for the cosmic 
ray spectrum, but here used directly from |54) . 

Following [53] , we assume that the production spectra 
for neutrons and cosmic rays from a single AGN are given 
by 



Q n (E n ,L p ) oc L p exp 

Q CT (E p ,L p ) oc L p exp 
where 



E n 



-E n 



E n 



E n 1 E b 1 (E n < E b ) 
E^ 2 (Eb < E n ) 

(1) 

E p ~ 1 E b ' 1 (E p < E b ) 
Ep 3 Eb (Eb < E p ) 

(2) 



• Q n and Q cr represent the neutron and cosmic ray 
spectrum respectively, as a function of the neutron 
and proton energies E n and E p respectively, 

• L p represents the proton luminosity of the source, 

• Ei, is the spectrum breaking energy which can vary 
from 10 7 GeV to 10 10 GeV for optically thick AGN 
sources, and finally, 

• B m ax is the cutoff energy beyond which the spectra 
fall off steeply. 

Using Eq. ([I]) the generic neutrino production spectrum 
from AGN's can be written as 



Q„ (£) « 83.3Q n (25£) 



(3) 



We now need to account for red-shifting in the ener- 
gies of the neutrinos propagating over cosmological dis- 
tances prior to arriving at the detector. It is convenient 
to describe the red-shifting in terms of the dimensionless 
red-shift parameter z, defined as 



A 
An 



= 1 + z. 



A and A being wavelengths of a propagating signal at 
detector and at source respectively. In terms of z the 



energy of a particle at source (Eq) and at the detector 
(E) can be related via 

En 

Thus, to account for red-shifting in the energy of the 
neutrinos we replace the source energy E in Eq. ^ by 
E(l + z). We now incorporate standard neutrino oscil- 
lations by multiplying the spectrum with the oscillation 
probabilities. The probability of a neutrino flavour v a 
oscillating to another vp is given by 

Pa->p = S a /3 

- 4 Tie (U* m U^U a] U; 3 ) sin 2 ' " » 

i>j 



AE 



However, as the distances involved are very large, oscilla- 
tion only provides a z-independent averaging effect over 
the three flavours. In all our calculations, unless other- 
wise mentioned 3 , the CP violating phase 5c p is kept 
and the 3a best-fit values of the mixing angles [55] are 
used, i.e., 



sm 2 (6» 12 ) = 0.321, sin 2 



= 0.47, sin 2 (6»i 3 ) = 0.003. 



The intensity at earth for an input spectrum 
Q [(1 + z) E, z] is given by 



1(E) 



(1 + z) 2 dU c dP gal 
47rd 2 dz dV c 



Q[(l + z)E,z]dz (5) 



with d,L and V c representing the luminosity distance and 
co-moving volume respectively. 

To obtain the maximal bound for the diffuse flux from 
the optically thick sources we start with Eb = 10 7 GeV 
in the input spectrum Q and carry out the above inte- 
gration using z m i n — 0.03 and z max = 6. The value of 
Eb is varied from 10 7 GeV to 10 10 GeV and the above 
integration is carried out for each case. The resulting 
I Ei (E) are then superposed to obtain the final bound. 
This is then normalised using the observed cosmic ray 
spectrum to give the upper bound of the diffuse flux for 
the three neutrino flavours at the detector. As may be 
expected, it leads to a result similar to that obtained in 
[54] . with, however, the results of standard oscillations 
incorporated. A related procedure is used for calculating 
the fluxes from optically-thin sources. We call this nor- 
malised upper bound of diffuse fluxes the MPR bound, 
and use this as the reference flux for all our calcula- 
tions. The MPR bound is a modification of the Waxman- 
Bahcall (W & B) bound [SB], where a uniform E~ 2 input 



In Sec. |Vl| we use the tribimaximal value of O23 = 45° that 
ensures perfect symmetry between t/ M and v T under standard 
oscillation. 
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spectrum of extragalactic cosmic rays was used to cal- 
culate the diffuse fluxes. This difference is noticeable in 
Fig. [2] where we have shown both these reference bounds. 
The resultant MPR bounds for both types of sources are 
shown in all the figures as unbroken gray lines. 



III. ROLE OF STANDARD NEUTRINO 
OSCILLATIONS IN RESTORING PARITY 
AMONG FLAVOUR SPECTRA 

As has been discussed in j39j , a flavour ratio of f : 2 : 
at source is reduced to the democratic f : f : f at the 
detector due to oscillations. In this section we study 
the effect standard oscillations have on spectral shapes 
and magnitudes of the neutrino fluxes, as they propa- 
gate from the source to the earth. We find that AGN's 
being very distant sources, standard neutrino oscillations 
play a very significant role in restoring equality among 
the three flavours in terms of not only magnitudes, but 
also the shapes of the diffuse flux at the detector. As a re- 
sult, though the flavour fluxes at some exceptional source 
might differ from each other in their spectral shapes sig- 
nificantly, neutrino oscillations during propagation en- 
sure that such differences are largely wiped out by the 
time they reach the detector. In addition to spectral 
shapes, oscillations also have the effect of bringing widely 
differing magnitudes close to each other (roughly within 
a factor of two). This implies that if an ultra-high-energy 
detector at the earth detects significant difference in mag- 
nitude and energy dependence among the AGN diffuse 
flux of the three flavours, it must be as a consequence of 
non-standard physics present in the oscillation probabil- 
ities during propagation. 

To demonstrate this, we assume two flavour spectra 
at source intentionally chosen to be widely differing and 
propagate them to the earth. We calculate the diffuse 
flux (to arbitrary normalisation) of the flavours arriving 
at the detector from all sources, assuming they give the 
same spectra for the two flavours at source. The result 
is shown in Fig. [T] We have checked that such a conclu- 
sion holds true in general as it does in this representa- 
tive case, and demonstrates that significant differences of 
shape or magnitude among the diffuse flux flavours if de- 
tected must be pointers to non-standard physics playing 
its role during propagation. 



IV. EFFECT OF NEUTRINO DECAY 

A. Introduction to neutrino decay 

Bounds on the life-times of neutrinos are obtained pri- 
marily from observations of solar (23] and atmospheric 
neutrinos. Observations from solar neutrinos lead to 

— > I CT 4 s/eV (6) 
m 2 



while, if the neutrino spectrum is normal, the data on at- 
mospheric neutrinos constrain the life-time of the heavi- 
est neutrino 

— > 1CT 10 s/eV. (7) 

In the following, we treat the lightest neutrino as stable 
in view of the fact that its decay would be kinematically 
forbidden, and consider the decay of the heavier neutri- 
nos to invisible daughters like sterile neutrinos, unparticle 
states, or Majorons. Neutrinos may decay via many pos- 
sible channels. Of these, radiative two-body decay modes 
are strongly constrained by photon appearance searches 
[57] to have very long lifetimes, as are three-body de- 
cays of the form v — > vvv which are constrained |58) by 
bounds on anomalous Zvv couplings 59J. Decay chan- 
nels of the form 

v { -> Vj + X (8) 
v X (9) 

where Vi represents a neutrino mass eigenstate and X 
represents a very light or massless invisible particle, e.g. a 
Majoron, are much more weakly constrained, however 
and are therefore the basis of our consideration in this 
section. When considering decays via the channel in 
Eq. Q we assume that the daughter neutrino produced is 
significantly reduced in energy and does not contribute 
to the diffuse flux in the energy range relevant for our 
purpose (1000 GeV to 10 11 GeV). A detailed study of 
the various possible scenarios for neutrino decay is made 
in [60]. 

Prior to proceeding, we would like to discuss cosmolog- 
ical observations of high precision which might be able to 
constrain models of decay via channels as in Eq. ^ in the 
future. These constraints are based on the determination 
of the neutrino mass scale as discussed in [63], or from 
the cosmic microwave background as discussed in [64] . 
Such observations would serve to push the lower bound 
of neutrino decay lifetimes by several orders of magnitude 
compared to those discussed here. However, these predic- 
tions are dependent upon the number of neutrinos that 
free-stream and assume couplings of similar nature and 
strength for all the species of the neutrino family. As dis- 
cussed in [65] and [66] these assumptions must await con- 
firmation and rely on future data. Hence, "fast" neutrino 
decay scenarios are not ruled out within the scope of cur- 
rent theory and experiment, though they are disfavoured. 
Further the decay of neutrinos via Eq. ([9| and in the 
cases where the decay, both via Eq. |8]) and Eq. ([9| hap- 
pen due to unparticle scenarios are not covered by such 
constraints and the purely phenomenological and general 
study of neutrino decay in the life-times discussed here 
would still be interesting and relevant for future neutrino 
detectors. 
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Spectrum for and v e at source Spectrum after propagation [z=2J 




E [GeV] 

FIG. 1: The even-ing out of possible spectral distortions present at source due to standard oscillations over large distances 
as seen for hypothetical spectra of two flavours (deep-red) and (green) from an AGN source at a redshift z = 2. 1(E) 
represents the flux spectrum for the two flavours. 



B. Effect of neutrino decay on the flavour fluxes 

A flux of neutrinos of mass m,, rest-frame lifetime Tj, 
energy E propagating over a distance L will undergo a 
depletion due to decay given (in natural units with c = 1) 
by a factor of 



<>Aj>; - //: t i = <<\j> ( -— x ^ 



where t is the time in the earth's (or observer's) frame 
and 7 = E/rrii is the Lorentz boost factor. This enters 
the oscillation probability and introduces a dependence 
on the lifetime and the energy that significantly alters the 
flavour spectrum. Including the decay factor, the proba- 
bility of a neutrino flavour v a oscillating into another vp 
becomes 



P aP (E) = J2\ u 



\U n 



|2 e -L/T s (£ 



Ka?P, (10) 



which modifies the flux at detector from a single source 
to 

<f>»«(E) = Y,<f^WPfii\ 2 \Uca\ 2 e- L frW.(iy 

We use the simplifying assumption T2/TO2 — T^/m^ = 
r/m for calculations involving the normal hierarchy (i.e. 
TOg — ml = Anij; > 0) and similarly, T\/m\ = t 2 /to 2 — 
r/m for those with inverted hierarchy (i.e. Am'j < 0), 
but our conclusions hold irrespective of this. The total 
flux decreases as per Eq. (11), which is expected for de- 
cays along the lines of Eq. ([9]) and, within the limitations 
of the assumption made in Sec. IV A also for Eq. p]). 

The assumption of complete decay leads to (energy 
independent) flux changes from the expected i/~ : lA : 
v% = 1 : 1 : 1 to significantly altered values depending 
on whether the neutrino mass hierarchy is normal or in- 
verted as discussed in [45]. From Fig. [2] we note that 
the range of energies covered by UHE AGN fluxes spans 
about six to seven orders of magnitude, from about 10 3 
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E [GeV] E [GeV] 

FIG. 2: Modification of MPR bound for incomplete decay with normal hierarchy (left) and inverted hierarchy (right), and 
life-time T2/rri2 = Tz/raz =0.1 s/eV. The z/ M and u e fluxes shown are from optically thick (in thick) and optically thin sources 
(thinner). Similarly the gray lines indicate the v e , v^, or v T undistorted flux modified only by neutrino oscillation, for both 
optically thick and thin sources, sensitivity thresholds and energy ranges of relevant experiments, viz., AMANDA and IceCube 
[61] , and ANITA [62] are indicated. 1(E) denotes the diffuse flux spectrum of flavours at earth, obtained as described in the 
text. 



GeV to 10 10 GeV. For the "no decay" case, the lowest 
energy neutrinos in this range should arrive relatively in- 
tact, i.e. L/E ~ t/itl ~ 10 4 sec/eV. In obtaining the last 
number we have assumed a generic neutrino mass of 0.05 
eV and an average L of 100 Mpc. On the other hand, if 
there is complete decay, only the highest energy neutrinos 
arrive intact, and one obtains i.e. L/E ~ r/m < 10~ 3 
sec/eV. Thus, a study of the relative spectral features and 
differences of flavour fluxes at earth allows us to study 
the unexplored range 10~ 3 < r/m < 10 4 via decays in- 
duced by lifetimes in this range (we have referred to this 
case as "incomplete decay" in what follows). 

To calculate the MPR-like bounds with neutrino decay 
we use the procedure of Sec. Ill] but replace the standard 



neutrino oscillation probability by P a p given in Eq. ( 10 ) 
with E replaced by E(l + z) to account for red-shifting. 
Since, unlike standard oscillations, P a p has an energy de- 
pendence that does not just average out, the diffuse flux 
obtained with decay effects differ considerably from the 
MPR bounds in shape as well as magnitude. Fig.[2]shows 
the effect for both normal and inverted hierarchies with 
a lifetime of Ti/mi = T^/m^ = 0.1 s/eV. We note that 
the effect of decay in altering the diffuse flux spectrum is 
especially strong in the case of inverted hierarchy. 

Fig. [3] shows how the diffuse flux spectral shapes 
change as the lifetimes of the two heavier mass- 
eigenstates are varied between 10 -3 s/eV and 1 s/eV. 
From the figure it is clear that this (10 -3 s/eV - 1 s/eV) 
is the range of life-times that can be probed by ultra- 
high-energy detectors looking for spectral distortions in 
the diffuse fluxes of the three flavours. For lifetimes above 
1 s/eV the spectral shapes start to converge and become 
completely indistinguishable beyond 10 4 s/eV while for 
those below 10~ 3 s/eV the shapes of the diffuse fluxes 
show no difference although their magnitudes are expect- 
edly very different. 



As is also the case for complete decays, the results are 
very different for the two possible hierarchies. This is be- 
cause the mass eigenstate m\ contains a large proportion 
of 24, whereas the state is, to a very large extent, 
just an equal mixture of and v T with a tiny admix- 
ture of v e . Therefore decay in the inverted hierarchy case 
would lead to a disappearance of the eigenstate with high 
content of v e and, hence, to its strong depletion against 
the other two flavours. In the normal hierarchy case, in 
comparison, the mass eigenstate with the high content 
of u e is also the lightest, and decay of the heavier states 
consequently leads to a depletion of and v T . Thus 
incomplete decay to the lowest mass eigenstate with a 
normal hierarchy (i.e. m,\) would lead to considerably 
more shower events than anticipated with an inverted 
hierarchy. 

While assessing the results presented here, it must be 
borne in mind that observation of a significant amount of 
V e from supernova SN1987A possibly imposes lower lim- 
its on decay lifetimes of the heavier neutrinos for the in- 
verted hierarchy scenario that are much higher than those 
considered here [S3 EH]- This observation, of a flux of V & 
roughly in keeping with standard predictions constrains 
its "lifetime" r/m > 10 5 , i.e., higher than what would 
give observable results with the methods described here. 
Despite the uncertainties involved with neutrino produc- 
tion from supernovae and the fact that the total signal 
from SN1987A was only a handful of events, the results 
for decay with inverted hierarchy must be judged keeping 
this in view. 



C. Modification of total UHE events due to decay 

The effect of decay as seen in the diffuse fluxes in Fig. [2] 
above must also translate to modifications in the shower 
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FIG. 3: Modification of MPR bound for incomplete decay with normal hierarchy (left) and inverted hierarchy (right), and 
life-times varying from r/m = 0.001 s/eV to 1.0 s/eV. The and v e fluxes shown are from optically thick sources. The gray 
lines indicate the z/ e , or v T undistorted flux modified only by neutrino oscillation. Similar effects are seen with fluxes from 
optically thin sources as well. 



and muon event rates observable at UHE detectors. In 
this section we demonstrate this by a sample calculation. 
We calculate the event-rates induced by the three flavours 
of high-energy cosmic neutrinos after decay using a sim- 
plified version of the procedure in Ref. |69j and compare 
it to those predicted by standard physics. 

Events at the IceCube will be classified primarily into 
showers and muon-tracks. Shower events are generated 
due to the charged current (CC) interactions of v e and 
v r below the energy of 1.6 PeV and neutral current (NC) 
interactions of all the three flavours. For energies greater 
than 1.6 PeV, CC interactions of the v r have their own 
characteristic signatures in the form of double-bangs, lol- 
lipops, earth-skimming events, etc. [70l[7T]. Muon-tracks 
are generated due to the induced CC events. 



To simplify Eq. ( 14 ) we use the approximation 



da 



at (y - (y)) 



(15) 



where (y) is the mean inelasticity parameter. Thus, we 
have 



poo 

Rnte = N A V dE u a NC (E u )T(E u ), (16) 

J E' th 

E' th is an effective threshold energy at which the curves 
defined by y — E t h/E v and y = (y) intersect. 



induced events 



v e induced events 

In the standard model v e interacts with nucleons via 
CC and NC interactions leading to electromagnetic and 
hadronic showers. 

In the CC events, the shower energy is equal to the 
initial neutrino energy E v , that is, the total energy of 
the two final state particles (an electron and a scattered 
quark). The event rate for v e N — > e~x, with \ being a 
final state quark, is given by 



Rate 



dy N A L*%° 
ay 



AT (E v ) 



N A V 



pOO 

/ dE v a C c{Eu)HE v ) 

JE th 



(12) 
(13) 



where 

• E v : the incident neutrino energy 

• Etu' detection threshold for shower events 

• y: the inelasticity parameter defined as y 



• A, L, V: the area, length and volume of the detector 
respectively 



• ?{E V ) 
GeV 



the flux spectrum of neutrinos in 



1 —2 —1 

cm s 



It is assumed that the electron range is short enough such 
that the effective volume of the detector is identical to the 
instrumental volume. Using standard tabulated values 
of the cross-section acc [ZH [73] it is straightforward to 
evaluate the integral in Eq. ( 13 ) to obtain the event rate. 



The event rate for anti-neutrino process V e N — » e + x is 
calculated similarly. 

For the NC events, the final state neutrino develops 
into missing energy, so that the rate is given by 



Rate = / dE v 

>E th 



dyN A L 



dffjyc 
dy 



AT (E v ) (14) 



The muon track event is calculated by 



OO pi- 

dE v N A \ 



dyR{E v {l-y),E., 



th ) 



x da ^S{E u )AF{E v ) 1 (17) 
dy 



where, 



R ^ v) = \ ln i^Vy 



(18) 



with a = 2.0 x 1(T 3 GeV cm -1 and b = 3.9 x 
10~ 6 GeV cnV 1 . S(E V ) represents the shadowing effect 
by the earth [721 175]. 



Approximating using Eq. (15 1 gives 



Rate= / dE v N A R{E v {l-{y)),E th ) 

x a C c{E v )S{E v )AT{E v ) (19) 

with E' th being determined similarly as for the v e induced 
events. 

Using the procedure described above, we calculate the 
total shower and muon-track detector events (for V + v) 
for the inverted hierarchy scenario with a life-time of 1.0 
s/eV depicted in Fig. [3] (top-right) and compare it to 
the events expected from standard physics. The results 
are tabulated in Table U where we show event rates for 
UHE detectors, like the IceCube, over a 10 year period 
integrated over solid angle. The difference between the 
ratio of muon-track to shower events due to standard 
oscillation and that after considering neutrino decay are 
shown in Fig. [4] 

The disappearance of a majority of shower events (due 
to the depletion of the v e flux compared to that of v^) at 
lower energies, and their reappearance at higher energies 
is a distinctive feature. It indicates the presence of new 
physics (like incomplete decay) as opposed to spectral 
distortions originating in the source, or the appearance 
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14 
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3 


3 
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imum affects the decay effected diffuse fluxes, while in 



TABLE I: Total shower and muon-track detector events (for 
V + u) over 10 years (rounded off to whole numbers), and 
integrated over solid angle for the inverted hierarchy scenario 
with a life-time of r/m = 1.0 s/eV depicted in Fig. [3] 



a 3 




FIG. 4: The ratio (R) of muon-track events to shower events 
with inverted hierarchy and life-time r/m — 1.0 s/eV as 
shown in Table [I] The ratio expected due to standard physics 
is shown in brown, while the modified ratio due to the effects 
of decay is shown in light red. At energies greater than 10 8 
GeV, R due to standard physics and that after considering 
decay become equal. 



of a new class of sources. In the latter case, a correspond- 
ing depletion and subsequent enhancement is expected in 
muon events. By contrast, in the case of incomplete de- 
cay the fluxes return to the democratic ratio at higher 
energies where the neutrinos do not decay. 



V. EFFECT OF NON-ZERO CP VIOLATING 
PHASE AND Bis VARIATION ON NEUTRINO 
DECAY 

As described in Sec. HVI the calculation for the effect of 
decay of heavier neutrinos on the diffuse flux spectrum 
was done keeping the CP violating phase Sep — and 
6*i3 at the 3cr best fit value which is close to zero. In 
this section we look at how our conclusions are affected 
if we change these parameters significantly. In Sec. |V A| 
we look at how changing 6*13 from to the CHOOZ max- 



Sec. V B we examine the consequences of a non-zero CP 
violating phase in the same context. 



A. Variation of 9i3 

Observations at CHOOZ [TJ] constrain the maximum 
value of 0i3 (90 % confidence level) such that 

sin 2 (26>?3 ax ) = 0.10. 

Therefore, we have for # 13 the following experimentally 
allowed range of values 

< 6>i3 < 9.1° 

We allow 0i3 to vary within this range and study its 
effect on the results of Sec. [W] The results are repre- 
sented in Fig. [5] ft is clear that the effect of varying 9 13 
is significant. However, given the strong difference in the 
diffuse flux spectra for inverted and normal hierarchies, 
variation of $13 over the entire range would not affect our 
qualitative conclusions in Sec. [TV] regarding differentiat- 
ing between the two. 



B. Non-zero CP violating phase. 

The CP violation phase in the three family neutrino 
mixing matrix is as yet not experimentally determined. 
Neutrino telescopes probing ultra-high energies might be 
able to improve upon our present knowledge of this pa- 
rameter (see [75] , for example) . Here we look at how the 
presence of a non-zero CP violating phase, Sep in the 



mixing matrix could affect results obtained in Sec. IV 

Sep enters the oscillation probability via the mixing 
matrix as the product sin (8 13 ) ■ exp (±iScp)- Therefore, 
a non-zero CP violating phase does not affect any of our 
calculations if $13 = and its effect is imperceptible even 
when the 3tr best-fit value of #13 is used as is the case in 
Sec. |IV| For the remainder of this section we keep 6> 13 at 
the CHOOZ maximum and vary the CPV phase from 
to 7r. Fig. [6] shows the result on the flavour for decay 
in the case of a normal hierarchy for diffuse flux from 
optically thick sources. In the same way Fig. [7] shows 
the effect of a non-zero CP violating phase on decay with 
both the normal and inverted hierarchy. The effect of CP 
violation is quite small on the diffuse flux with inverted 
hierarchy as compared to that with normal hierarchy. 

To summarise, it is clear from the discussion in Sec. |IV| 
and [V] that future neutrino detectors capable of distin- 
guishing between flavours should be able to probe and 
potentially provide stronger bounds on decay lifetimes of 
heavier neutrinos. If the neutrinos decay with a lifetime 
within the ranges discussed here, then they would also be 
able to distinguish between the two hierarchies due to the 
strongly different diffuse flux spectra the two hierarchies 
lead to for the flavours v e and Vn, notwithstanding the 
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FIG. 5: Effect of variation of #13 over the complete range on decay plots obtained in Fig. [2] using optically thick sources. 
The shaded regions indicate the area spanned by the diffuse flux spectra as #13 varies from to the CHOOZ maximum, while 
the thick lines represent the spectra obtained with the 3a best-fit value of #13. 




E |GeV] 



FIG. 6: Effect of CP violation on the diffuse flux of the i/ M 
flavour obtained by considering decay with normal hierarchy 
and life-time of r/m = 0.1 s/eV. The variation in the flux as 
the CP violating phase is varied between - n is shown as 
the shaded region. 



effect of a non-zero CP violating phase or the uncertainty 
over the value of #13. 



VI. EFFECT OF LORENTZ SYMMETRY 
VIOLATION 

Low energy phenomenology can be affected by Lorentz 
symmetry violating effects originating at very high ener- 
gies. Typically such effects originate at energies close to 
the Planck scale. They may appear in certain theories 
which are low energy limits of string theory [251 176) , or 
could possibly signal the breakdown of the CPT theo- 
rem |77j . Additionally, if quantum gravity demands a 
fundamental length scale, leading to a breakdown of spe- 
cial relativity, or loop quantum gravity [76HF3] leads to 



discrete space-time, one expects tiny LV effects to per- 
colate to lower energies. UHE neutrinos, with their high 
energies and long oscillation baselines present a unique 
opportunity for testing these theories. Their effects in 
the context of flavour flux ratios have been discussed in 
|50) . They may arise, for example, due to a vector or 
tensor field forming a condensate and getting a vacuum 
expectation value, thereafter behaving like a background 
field. The effective contribution of such background fields 
can then be handled in the low energy theory using stan- 
dard model extensions [23]. It has been shown [77] that 
although CPT symmetry violation implies Lorentz vio- 
lation, Lorentz violation does not necessarily require or 
imply the violation of CPT symmetry. In this section we 
focus on the modification of the propagation of neutri- 
nos due to Lorentz symmetry violating effects along the 
lines discussed in Ref . [24] . Since the effects of Lorentz- 
violation and CPT violation are understandably tiny at 
low energies, it is difficult to explore their phenomeno- 
logical signatures using low energy probes, in colliders 
for example. Since they originate in extremely energetic 
cosmological accelerators and propagate over cosmic dis- 
tances, ultra-high energy neutrinos provide the perfect 
laboratory for constraining and, possibly, determining 
Lorentz- violating parameters. 



A. Modification of neutrino transition probabilities 
due to LV effects 

As an example, we will study, for the simplification 
that it provides, a two-flavour scenario with massive neu- 
trinos and consider the modification of the transition 
probability from one flavour to the other by Lorentz- 
violation due to an effective standard model extension. 
Our focus is on LV from off-diagonal terms in the effective 
hamiltonian describing the propagation of the neutrinos 
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FIG. 7: Effect of CP violation on fluxes affected by decay for both normal and inverted hierarchies. The shaded regions 
represent the span of the flux bounds when the CP violating phase is varied from to ir, keeping the #13 at the CHOOZ 
maximum. 



We consider an effective Hamiltonian describing neu- 
trino propagation 



1 



H% HP I <W + 2T^T + 2 (0%)] at 



(20) 



where m is related to the neutrino mass and a is a real 
CPT and Lorentz violating parameter. In the two neu- 
trino mass basis this gives 



with 2 3 = 45° is a unique signature of Lorentz- violation. 
This would lead to a significant decrease in the signature 
v T events at high energies, like "double-bang" , "lollipop" 
and "earth-skimming" events as compared to muon-track 
events. Differences in shape between the two flavours 
can be seen for a < 10~ 30 GeV. We have used the case 
where a is independent of energy, however if the param- 
eter a oc E n the results would be qualitatively similar 
to that obtained here but involve significantly different 
ranges of values for the parameter as expected. 



'off 



(21) 



With the mixing angle between the two flavours #23 = 
7r/4, this modifies the probability of transition from one 
flavour to another during propagation to 



where u> = 



Am" 
4E 



and £1 



1 



a 

fi2 



^cos(20L)j 



(22) 



y/to 2 + a 2 . 



B. Effect of Lorentz violation on neutrino flavour 
fluxes 

To calculate the diffuse fluxes of the two neutrino 
flavours we use Eq. ( 22 1 instead of the standard os- 



cillation probability and integrate over the red-shift z. 
The probability above contributes a z dependent term 
through its dependence on energy. Further the cos (2f2£) 
term averages out and consequently does not contribute. 

The results of including Lorentz violation in the prop- 
agation phenomenology of neutrinos are shown in Fig. [8] 
It is clear from these plots that the strong departure of 
diffuse spectral shapes of and v r from the symme- 
try expected under standard oscillation phenomenology 



C. Detectability of Lorentz- violation 

Unlike in neutrino decay, the effect of Lorentz violation 
is seen in the deviation of the flux spectra of both the 
and, more strikingly, the v T flavour, from the standard 
fluxes toward the higher end of the spectrum. This makes 
it especially interesting for probe by detectors, such as 
ANITA and the Pierre Auger Observatory [SU [55] hav- 
ing sensitivity to v T in the energy range 10 s — 10 11 GeV. 
While Auger can separate out the v T events, ANITA de- 
tects the sum of all three flavours. As is clear from the 
experimental thresholds shown in Fig.JHJ should even tiny 
Lorentz- violation effects exist, both these experiments 
will, in principle, be able to detect it via lack of char- 
acteristic r events expected at these energies from stan- 
dard physics. As they collect more data in the future, 
expectedly bringing the corresponding thresholds down, 
the ability of such experiments to detect tiny LV effects 
will be gradually enhanced. 



VII. PSEUDO-DIRAC NEUTRINOS 

Masses for neutrinos can be generated by extending 
the Standard model to include right-handed sterile neu- 
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FIG. 8: Effect of Lorentz violation on the — v T diffuse flux with various values of the lorentz violating parameter a (in 
GeV). Clockwise from top- left (i) a = 0, (ii) a = 10~ 30 , (iii) a = 10 -28 , (iv) a — 10~ 26 . The plots show how an increase in 
the LV parameter results in depletion of the v T flux at progressively lower energies. For the Auger experiment, sensitivities for 
v T detection using the most pessimistic systematics (top line) and the most optimistic systematics (bottom line) are indicated 



trinos to the particle spectrum. The generic mass term 
for neutrinos becomes 

£ = --tOM$ + /i,c., (23) 

where considering 3 right-handed neutrinos in the spec- 
trum 

* = {vcL, v^l, v tL , {v 1R ) C , (v 2 r) C , (^3fl) C ) , 

and v G = Cv 1 ', C being the charge conjugation operator. 
The mass matrix M is of the form 

M = ( mL m ?) , (24) 
\m D m R J 

and for = mp , = reduces to neutrino states with 
Dirac mass. In this case the six neutrinos decompose into 
three active-sterile pairs of neutrinos degenerate in mass 
with maximal mixing angle 9 — 7r/4 for each pair. Due to 



the mass degeneracy within the neutrinos in such a pair, 
an active neutrino cannot oscillate into a sterile neutrino 
from the same pair. 

Instead, neutrinos may be pseudo-Dirac states |48j 
where wil and mn are tiny but non-zero, i.e. ttil, tur <C 
mjj. This lifts the degeneracy in mass within an active- 
sterile pair, and gives a mixing angle 9 7r/4 between its 
members. The result of the lifting of this degeneracy is 
to enable oscillation among species that was not possible 
in the pure Dirac neutrino case. 

The presence of non-zero rax, , m r changes the proba- 
bility of transition of one active state to another during 
propagation. The expression for the probability for neu- 
trinos propagating over cosmological distances (after var- 
ious phase factors involving terms like Atoq/L average 
out) is g8] 

P a0 = £ I U aj | 2 | U Pi | 2 cos 2 — >— , (25) 
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FIG. 9: Effect of pseudo-Dirac (PD) neutrinos on the 

v a diffuse flux with Am 2 = 10" 14 eV 2 . 



where Am 2 = (m^~) — (m^ ) 2 is the mass squared dif- 
ference between the active and sterile states in the j th 
pair. 

There has been a recent study [55] that explores the 
pseudo-Dirac scenario at neutrino telescopes using the 
ratio of shower to muon-track events. Here, we look at 
distortion of spectral shape from the standard diffuse flux 
due to the modification of the oscillation probability to 



Eq. (25 1. We use Eq. (25 1 instead of the standard oscilla- 



tion probability, otherwise following the same procedure 
used to derive the standard MPR flux (the base flux in 
our plots). The results are shown in Fig. [9] which shows a 
decrease in the affected flux at lower energies and rise at 
the higher end of the spectrum to merge with the stan- 
dard flux. However, the decrease is only to about half 
the base flux, and the rise at higher energies is not steep. 
Therefore, it would be very difficult to detect such an 
effect in future detector experiments. 



VIII. EFFECT OF DECOHERENCE DURING 
NEUTRINO PROPAGATION 



Quantum decoherence arises at the Planck scale in the- 
ories where CPT invariance is broken independently of 
Lorentz symmetry due to loss of unitarity and serves to 
modify the time evolution of the density matrix [49j [50] . 
Though not expected in a majority of string theories, a 
certain class of string theories called noncritical string 
theories may allow for decoherence. 

In the context of neutrino oscillation, decoherence 
serves to modify the transition probabilities among the 
three flavours. While a general treatment discussing how 
this happens for the three family case is complicated, we 
work under the simplifying conditions assumed in |501 see 



Sec IV. B] to arrive at the transition probability 



P [v p v q ] 



1 1 



er 2SL [3 (U 2 pl - U 2 2 ) (U 2 ql - U 2 q2 ) 



3 ' 6 

+ (C/pi + U 2 2 - 2U 2 3 ) (U 2 ql + U 2 2 - 2U 2 3 )] , 

(26) 



where S is the only decoherence parameter. This leads to 
a flavour composition at the detector given by 



Rv, 



4>, 



TOT 



$T 



R v =P\Ve 



$T 



^T 



(27a) 



P [V T -> V,x\ 



<J> TOT 



(27b) 



PW 



PWr 



TOT 



TJ ^T 



(27c) 



where <I> e /$ T oT, etc. are flux composition ratios at source. 



We use the flavour ratios given by Eq. (27) to calcu- 



late the diffuse flux spectra of each flavour arriving at 
the detector. The effect of decoherence is to bring the 
flavour fluxes close to the ratio 1:1:1. If we use the 
standard flux from AGN's (1 : 2 : at source) then stan- 
dard neutrino oscillation already brings the ratio to the 
above value as discussed in Sec. lIIII and this makes it diffi- 
cult to distinguish between the effects of decoherence and 
standard oscillation. However, if we have detection capa- 
bilities that can distinguish between neutrinos and anti- 
neutrinos, it might be worth investigating decoherence 
using the differences in flavour spectral shapes. As dis- 
cussed earlier pion decays in the source via 7r + — > 
and subsequently, /i + — > e + V fl v e contribute to a flavour 
spectral ratio of 1 : 1 : for v and 0:1:0 for V. Due 
to standard oscillation these flavour ratios are reduced to 
0.78 : 0.61 : 0.61 and 0.22 : 0.39 : 0.39 at the detector re- 
spectively. Since the effect of decoherence is to reduce the 
flavour ratios to 1 : 1 : 1 irrespective of ratios at source, 
the transition from the flux due to dominance of stan- 
dard oscillation to that due to dominance of decoherence 
might happen within the energy range relevant for our 
purposes, for a certain range of values of the decoherence 
parameter. However the effect is almost invisible even if 
v and v fluxes are used as probes, the reason being that 
the fluxes ratios at detector due to standard oscillation 
for both (i.e., 0.78 : 0.61 : 0.61 and 0.22 : 0.39 : 0.39 
respectively) are already quite close to the 1:1:1 that 
decoherence would result in. Effective probe for decoher- 
ence are high energy neutrinos from neutron decay, for 
instance, which gives a flux ratio of 1 : : at source |69j . 
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FIG. 10: Effect of decoherence on the diffuse flux with the parameter 5 = aE 2 and a = 10 -40 GeV -1 . A base flux 
composition of : 1 : corresponding to V (left) and 1:1:0 corresponding to v (right) from pion decay is used for the 
calculation. It is clear from the figure that (anti-)neutrinos from pion decay are not useful probes for decoherence. 



and not neutrinos from pion decay. The results for V and 
v with a particular choice of the decoherence parameter is 
shown in Fig. |10| For our calculation, we have chosen the 
parameter S oc E 2 which is expected within the context 
of string theories 4 . Upper limits on such a parameter are 
got from the Super-Kamiokande as ~ 1CP 10 GeV. 



IX. CONCLUSIONS 

In this article we have discussed the effects of several 
exotic, non-standard physics on the diffuse fluxes of the 
three neutrino flavours, using neutrino fluxes from AGN's 
as an example. We have assumed a standard neutrino 
flux at source with the flavour ratio thereof being 1:2:0 
and shown that due to standard oscillations in vacuum 
during the propagation of these neutrinos across cosmo- 
logical distances the fluxes are evened out to the demo- 
cratic value of 1 : 1 : 1, and that even for non-standard 
fluxes at source the fluxes at the detector are still close 
to each other in magnitude and their spectral shapes are 
very similar. 

Non-standard physics serves to destroy this equality 
among the three flavours and this serves as a potential 
probe for the underlying nature of the physics involved. 
To demonstrate this we first looked at how the decay of 
the heavier of the neutrinos affects the standard MPR dif- 
fuse flux bounds in the case of both normal and inverted 
hierarchies. We found that decay life-times of magnitudes 
several orders above those currently understood from ex- 
periment induce detectable changes in spectral shapes of 
the three diffuse fluxes, both against the standard flux, 
and among each other. Since the effects are strikingly dif- 



4 The choice of <5 oc E 2 also violates Lorentz symmetry which 
introduces weaker secondary effects not taken into account here. 



ferent for the two hierarchies, it would also be possible to 
search for the hierarchy in case the heavier neutrinos do 
decay with life-times in the range 10 -3 s/eV - 10 4 s/eV, 
as discussed here. We have also shown that the effects 
remain significant despite variation on the unknown pa- 
rameters 6*i3 and Sep and probing neutrino decay within 
the life-times explored here should be possible despite our 
limited knowledge about these parameters. 

Tiny effects of Lorentz symmetry violation in the low 
energy theory arising due to the effect of some Planck 
scale physics can also be probed using ultra-high energy 
neutrinos. Taking the simplifying case of two neutrino 
flavours and v T we have described the effect of Lorentz 
violating parameters on transition probabilities between 
them during propagation and inferred that it leads to a 
strong decrease in the v T flux as compared to the v„ flux. 
This breaks the — v T symmetry that is a feature of all 
standard model and most beyond standard-model sce- 
narios, and thus provides us with a distinctive signature 
for LV. It translates to a corresponding decrease in the 
signature v T events at high energies. While a simplifying 
case of two flavours and involving just the one Lorcntz- 
violating parameter was dealt with here, the conclusions 
are true more generally. Detection of a sharp decrease 
in r events in future detectors like Auger and ANITA 
will be an indicator of the extent of Lorentz violation in 
low energies. Conversely, the failure to detect such a dip 
could be used to put bounds on the LV parameters. 

Further, we have discussed the effect of decoherence 
and the existence of pseudo-Dirac neutrino states on the 
diffuse fluxes of the three flavours. While not as strik- 
ing as the effects of neutrino decay or LV, the existence 
of pseudo-Dirac states affects distortions in the spectral 
shape of the standard flux at the lower end of the spec- 
trum. On the contrary, decoherence shows almost no 
distortion on the fluxes. A probe of decoherence requires 
that we distinguish between neutrinos and anti-neutrinos 
since, irrespective of the flux ratio at source, it tries to 
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bring the flux ratio to 1 : 1 : 1 at the detector, same 
as what standard oscillation does to the standard flux of 
1:2:0. Even so, the effect of decoherence, seen at higher 
energies, is not significant and cannot, in all probability, 
be experimentally distinguished. 

It is clear that future ultra-high energy neutrino detec- 
tors with strong flavour detection capabilities and excel- 
lent energy resolution will allow us to probe the validity 
of non-standard physical phenomena over large ranges 
of the involved parameters. While differences in spectra 
among the flavours arise due to the selectivity of non- 
standard physics with regard to the three families, strong 
distortion of spectral shape of the fluxes as compared to 
the standard flux expected at the detector arises due to 
the non-trivial energy dependence of transition probabil- 



ities in new physics. To detect or, potentially, constrain 
new physics it is necessary to carry out experiments that 
combine searches of both kinds. While understandably 
challenging, it will certainly be worthwhile carrying out 
detection experiments along these lines given the funda- 
mental nature of physics that will be brought under the 
scanner. 
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